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We formulate higher derivative gravity with Lagrange multiplier constraint and scalar projectors. 
Its gauge-fixed formulation as well as vector fields formulation is developed and corresponding 
spontaneous Lorentz symmetry breaking is investigated. We show that the only propagating mode 
is higher derivative graviton while scalar and vector modes do not propagate. Despite to higher 
derivatives structure of the action, its first FRW equation is the first order differential equation 
which admits the inflationary universe solution. 
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I. INTRODUCTION 

The study of the evolution of early-time and late-time universe indicates that General Relativity is not complete 
theory of gravitational interaction. At best, it may be the effective theory to describe the classical gravity at in- 
termediate energies. Indeed, there are quite strong proposals on the unified description of the early-time inflation 
and late-time acceleration in terms of modified gravity (for review, see [1]). It is desirable to extend such modified 
gravity till Planck scale where it should be renormalizable or somehow consistent (finite?) one at UV. Otherwise, 
such proposals which pass local/cosmological tests (for recent review, see remain to be phenomenological ones. 

It is quite well-known that higher derivative gravity may be multiplicatively renormalizable (for a review, see 
Q). However, the unitarity issue in such approach is so far the open problem. Recently, very interesting attempt 
to formulation of power-counting renormalizable gravity has been presented in terms of gravity [3] which also uses 
effectively higher derivative propagator. Unfortunately, such approach leads to explicit breaking of Lorentz invariance 
from the very beginning. Nevertheless, such formulation may be extended to covariant theory Q where Lorentz 
invariance is broken spontaneously, in the same way as celebrated gauge invariance. The construction of higher- 
derivative power-counting model with consistent graviton spectrum in this direction may suggest the new perspectives 
towards to multiplicative renormalizability as well as resolution of some problems which appear in Lorentz non- 
invariant gravities. 

In the present work we propose covariant Lagrange multiplier constrained higher derivative gravity with scalar 
projectors. Its gauge-fixed formulation is developed. It is demonstrated that only higher derivative graviton degree 
of freedom is propagating at tree level while scalar and vector degrees of freedom do not propagate. The theory turns 
out to be power-counting (super-)renormalizable. The spontaneous Lorentz symmetry breaking is investigated. The 
equivalent representation in terms of vectors is also developed. Finally, FRW cosmology is studied. It is shown that 
theory looks very similar to gravity. The possibility of eternal inflation is demonstrated. 



II. THE MODEL OF POWER-COUNTING RENORMALIZABLE, COVARIANT HIGHER 

DERIVATIVE GRAVITY 

It is expected that power-counting renormalizable covariant gravity should be higher derivative theory, for instance, 
of the sort recently proposed in refs. [5;]. In the present section we propose new formulation of such theory and derive 
its gauge-flxed formulation and propagator. It will be shown that it may have very good UV behavior in gauge-fixed 
formulation. Different versions of Lagrange multiplier gravity were discussed in refs. jlj]. 

In general, the unitarity is broken in the higher derivative theories. In order to guarantee the unitarity, the models, 
where Lorentz symmetry and/or the full general covariance is explicitly broken, were proposed by Hofava 14|. Although 
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the models are very interesting, due to the lacking of the Lorentz symmetry and/or the full general covariance, there 
appears an extra propagating and problematic scalar mode (see Q, for example). In the models proposed below, there 
is the Lorentz symmetry and/or the full general covariance in the actions but the symmetry and/or the covariance 
is broken spontaneously. As a result, we obtain the models where the UV behavior of the graviton propagator is 
improved but any extra mode, like scalar mode, docs not appear. 

We now start with the action including the Lagrange multiplier field A and the scalar field <j): 



Shag 

which gives a constraint 



d''xV^x(^^d^cl)d^<P + Uo^ , (1) 

^d^<l>d^cl> + Uo = 0, (2) 

that is, the vector is time-like. Therefore the Lorentz symmetry and/or the full general covariance is broken 

spontaneously. The detailed discussion about the spontaneous breakdown of the symmetry and/or the covariance will 
be discussed in detail in the next section. For the spontaneous breakdown, Uq needs not to be a constant but some 
non- vanishing and positive function of 4>- Just for the simplicity, we consider the case that Uq is a constant. At least 
locally, one can choose the direction of time to be parallel to {d^(f))- Then Eq. ^ has the following form: 

Therefore the spatial region becomes a hypersurface where is a constant since the hypersurface is orthogonal to the 
vector (9^0). 

We are interesting in the fluctuations over the fiat background: 



Note that ([3]) gives 

Then one can define a projection operator 



ff/xi. = fitly + h^iu . (4) 
(j, = A/2[^t . (5) 



and it follows 



P/P/i?^. = \ {h",^,, + y,^,, - d^d^h,,, - d,d, (V)) , (7) 
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and therefore 



— ''P/^p^VV^^V.^ = ihu,jt + htJ,^t - h,,,tt) , (8) 
a'^^a^V^V, + 2[/oV''Vp = 2Uodkd'', (9) 

P^ "P, " (Pm^ - ^^p^VV^V.^) ^ i {h,,/ + - h^^/ - d.d, (/j/)) . (10) 

Note that Pg^ = 0. Then we can propose the action of the power-counting renormalizable, covariant higher derivative 
gravity with scalar projector as^ 



'2n+2 



(11) 



We have merely found the expressions of the actions JTTJ and 1121 1 so that the propagator of the graviton could be given in 14711 . We 
have not found any deep physical principle to choose the actions in the forms of Ijllll and 11211 but we believe that there might exist 
some hidden symmetry to construct such theories. 



3 



for z = 2n + 2 model (n = 0, 1, 2, ■ 



and 









1 dt^xy 




.2^-"l 



(12) 



for z = 2n + 3 model (n = 0, 1, 2, • • • ). Here the quantity z is introduced to express the anisotropy between the time 
coordinate and spacial coordinates in (Q). 

One can confirm that the actions admit a flat space vacuum solution. Indeed, fleld equations are: 



A. 



(13) 



Here Cj^'^hcr ^Qj-^es from the higher derivative term (the second term) in the actions (|lll) and (IT2|) . Assuming the 
flat vacuum solution, Eq. ([2]) given by the variation over A has a form ([3]). Since for the flat space solution all the 
curvatures and V^Vi,0 vanish, Eq. ([T^ reduces to 

= A9^(/.a,</>, (14) 

whose solution is A = since 8^4) ^ due to the constraint equation ^ (if the coordinate system is chosen properly, 
we have dt(t> = -\/2[/o and di(l) = 0). Hence, the actions pT|) and (|12p admit the flat space vacuum solution with A = 0. 
Let us investigate the perturbation from the flat background ^ with A = in more detail. 

First, by using the diffeomorphism invariance with respect to time coordinate, we choose Eq. ([5]) as a (unitary) 
gauge condition. Then the actions (|lll) and (IT^ have the following form: 



52n+2 



-5? + dkd'') 



J i jk , J j ik J ij k 



S2n+3 — 



(15) 



{-2hu {s^'dkd'^ - d'd^) + 2hu {s'^dkO'^ - d'd^) htj 



-2^--'aU^"+' {(9.9^)" {h,,/ + h,^/ - h,^/ - 9.9, 



^dl + 9fc9'=) 



X { (9,9'^)"+' [K + V - h^\k - d'Q' {K)) } + C^oA/i«' 



(16) 



Here, only the terms quadratic with respect to the perturbation are kept. Note that there remains the diffeomorphism 
invariance with respect to the spatial coordinates. It will be fixed later in (|25p. An important thing is that there does 
not appear hu in the higher derivative term with a coefficient a. The constraint equation ([2]) shows 



htt = . 



(17) 



The idea proposed in ref. for quantum gravity is to modify tfie ultraviolet behavior of the graviton propagator in Lorentz non-invariant 
way as 1/ |fcp^, where fc is the spatial momenta and z could be 2, 3 or larger integers. They are defined by the scaling properties of 
space-time coordinates (a;, t) as tc — > hx and t — >■ h^t. When 2: = 3, the theory seems to be UV renormalizable. 
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The variation of hu can be solved with respect to A: 



for the action pTj) and 

A = ('5^^3,9'^ - 5^9^") /i,, + 2'-aU^- (9.9'^)'"+' a^c)-'" (/i,,/ + h,^/ - h,^/ - 9.9, • (19) 

for the action (jl2p . The hnearized equations given by the variation over cf) are: 

= 9i {a + 22«-iat/o'"-i (9fc9^f " 9^9-' (/ifc,/ + /i,,- >■ - h,^/ - 9,9, (/i^')) } , (20) 
for the action (fTTj) and 

= 9t {a + 22"aC/o2" (9^9'=)'"+' 9^9^ (/i,,/' + /i^^./ - h^^/ - 9,9, (/./))} . (21) 



for the action (|T2|). 

We now decompose /it,, which corresponds to the fluctuation of the shift function TVj, as follows 

hu = 9,s + V, , d'v, = . (22) 

Here s is the spatial scalar. We further write the linearized diffeomorphism invariance transformations with respect 
to the spatial coordinates as follows 

Sx' = d'u + w' , d,w' = . (23) 
Then under the diffeomorphism, s and Vi in p2p are transformed as 

Ss = dtu , 5vi = dtWi , (24) 
and therefore one can choose the gauge condition s = v"^ ~ 0, that is, 

hu^Q. (25) 

Furthermore the variation of hu gives 

dt {-25^^d' + + hjk = , (26) 

which is identical with that in the usual Einstein gravity since the higher derivative terms in the actions (|15p and (jl6p 
do not contain hu. We decompose hij as 

/ly = S,,A + 9, A: + 9,B, + a, + (^9,9, - y^J^k^''^ E , (27) 

with 

d'Bi = , 9*Cy = 9^C,,- = , = , (28) 
Then substituting (|27|) into ([26]). we obtain 

= 9* (^-49,A + 29fc9'=B, + ^dAd'^E^ . (29) 

By multiplying with 9% one gets 

9*9,9* (-AA + ^9fe9'=i;') = , (30) 



which shows 



A=i9fe9'=i?, (31) 



which also indicates that Bi = under the boundary condition that Bi should vanish at spatial infinity. 
Using ([27]) with Eqs. ^ and ^ have the following forms: 



A 



1 



1 



-A + -djd'E ) - 22"aC/2"-i {dkd 



under the boundary condition that A and E should vanish at spatial infinity. Then the equation given by substituting 
HI) into (gni) gives 

(32) 

(33) 
(34) 

(35) 
(36) 

(37) 



= a* |a + 2^"au^^-' {dkd'^y"^' (^-A + ^d.d^E^ 



and (HI]) and dm 



-dkd'^ i -A 



2k^Uq 

Combining (1311) with the above equations, we find 



djd^E] -22"+iaC/2" {dkd' 



-A + -djd^E 



-A + Id.d'E 



X^B^^O, 



and therefore the scalar modes A and the vector mode Bi do not propagate. 

After the gauge fixing and using ([571) etc., the actions ([TT]) and (IT^ have the following forms: 
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5. 



2ri+2 



1 

8^ 



{a, {-df + dkd'') c^^} - 2'-'au^- {{dkd^y^' a,} {{dkd'^y''' c^^} 



_L |_6A (-92 + dkd'') A + ^dkd^E (-92 + dkd") dkd^E 

AAdkd'A + ^A {dkd^^E - ^dkd^E {dkd'^f E^ 

22"-2a[/2" 1(9^9''^)" (^~d,djA ~ S.^dkd^A + ^d,d,dkd''E + ^5,, {dkd") 
I {dkd'')" (^-d'd^A - 5'^ dkd'' A + ^d'd^dkd''E + {dkd''f E 



E 



<S'2n+3 = / d'^X 



+^ |-6A (-92 + dkd") A + ^5^9'= £; (-92 + d^Qk^ d^d"E 
+4Adkd"A + ^A {dkd")^E - ldkd"E {dkd")^ E^ 

-22"-iat/2"+^ I {dkd")" (^-d,d,A - S,,dkd"A + ^d,djdkd"E + ^S,, {dkd")^ E 
d'd^A ~ 6'^ dkd" A + ^d'd^dkd"E + {dkd")^ E^ | 



{dkd") 

Then by the variation of A, we obtain 

= +^ |-12 {~df + dkd") A + 8dkd"A + ^ {dkd")^ E^ 

r,2n-l„,TT2n 



(38) 



(39) 



at/o^" {-d.,dj - 5,jdkd") I {dkd"f" (^d'd^A - 6'^ dkd" A + ^d'd^dkd"E + ^5'^ {dkd"f E^ | , 



(40) 



^ It is interesting that the structure of the actions II38I I and II39I I reminds the one of (7(1) invariant F(R) Hoi'ava-Lifshitz gravity [^. 
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for the action (l38l) and 



= +^ |-12 {-dl + dud'') A + md^'A + ^ (9^9'^)' £;| 



(41) 



for the action pQ]) . On the other hand, by the variation over E, one gets 



= a.a''^ 



i2rt-l 



+^^— (-9,9, - 5,,dkd'') I [dkd'^f' ( ^d'd^A - ,5^^ 9^9''^^ + ^9^9^9fc9'=i; + ^(5^^ (9^9*^)' £; 



(42) 



for the action 



and 



= 9fc9' 

)2n 



^ I ^ {-dl + 9fe9^) 9^9^=^ + ^dud'^A + ^ (9,9'^)' i?| 
+^aC^o"^' (-5»5j - %afc9^) I (9fc9'=)'"+' (^-9'9^A - <5'-'9fc9'=A + ^d'd^dud'^E + ^5'^ (9fe9'=)' 



for the action By using both of Eqs. (gOl), (HB), (HH), and (gS]) give the same expression: 



= dtA. 



(43) 



(44) 



Therefore A and E only depend on the spacial coordinate and therefore they do not propagate. Then we have shown 
that all the scalar modes A, hu, s, A, and E and all the vector modes Vi and Bi in Eqs. ((^ and (l?7l) do not 
propagate and the only propagating mode is massless graviton corresponding to the tensor mode Cij, which should 
be distinguished from Hofava quantum gravity Q where Lorentz invariance is explicitly broken. Then the actions 
([Tl]) and (fT2|) and therefore (|38| and (|39|) reduce to the simple forms: 



S2n+2 — 



2n+3 



J d^x {a, {-dl + dud^) C'^} - 22"-2aC/2» {(9,9'=)"+' C,, } { (9,9'=)"+' C^=] 
^ {a, (-9,2 + dud^) C'O - 22-VC/2-+1 {(9,9'=)"+' a,} {(9,9'=)"+' C^=] 



(45) 
(46) 



In the original models due to the traceless and transverse conditions for dj in ([28l) . the higher order terms do 
not contribute to the propagator of the graviton and therefore the UV behavior has not been really improved. 
Thus, the propagator has the following form in the momentum space: 



Stk - 



Pipk 

p2 



P J \ P J 

(p2 _ 22"-iaK2f^2n+ip2(2„+3))-^ ^ z = 2n + 3 case 



p2 / 



p. 



p' J 



z = 2n + 2 case 



(47) 



Here = Si=i {P^Y ^^'-^ 



p2 If a > 0, there appears the tachyonic pole when 



z = 2n + 2 case 



1 = 22"aK2c/2„p4„+2 ^ 
1 = 22"-iaK2C/2"+V*("+'^ , z = 2n + 3 case 



(48) 



with p'^ — and therefore at least the flat vacuum becomes unstable. On the other hand, there exist a stable flat 
vacuum when a < 0. 

In the present model, there is no propagating vector or scalar mode at least on the tree level. The change of the 
tensor structure of the propagator in (j47|) means that the vector or scalar mode could appear, that is, the vector or 
scalar mode must correspond to a composite state, which usually does not appear at any perturbative level. Therefore, 
it is expected the tensor structure should not be changed by the quantum corrections. 
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In the ultraviolet region, where k is large, the propagator behaves as 1/ \kf for z = 2 (n = 0) case in (|45|) and 
therefore the ultraviolet behavior is rendered. For z = 3 (n = 0) case in the propagator behaves as 1/ |fc| and 
therefore the model becomes power-counting renormalizable. For z = 2n + 2 (n > 1) case in (I45|) or z = 2n + 3 (n > 1) 
case in (|46p . the model becomes power-counting super-renormalizable. The dispersion relation of the graviton is then 
given by 

w = cofc^ , (49) 

in the high energy region. Here Cg is a constant, to is the angular frequency corresponding to the energy and k is the 
wave number corresponding to momentum. If cq < 0, the dispersion relation becomes inconsistent and therefore cq 
should be positive. 

Let us discuss the generality of the expression (|10p , which appears in the actions pT|) and . We now require for 
this expression, when a flat background is chosen 

1. The expression does not vanish. This condition is trivial but this often occurs due to some identity and the 
condition without torsion V fj,gup- 



2. The expression is given by the second rank symmetric tensor as in (jlOp. 

3. Each term contains the second derivative of perturbed metric h^^, like dpdah^j^^. This condition plays a role 
of the dimensionality condition. Since can be regarded as a constant vector, the third derivative of like 
V^iVy Vp(^ contains the second derivative of the perturbed metric like Opdah^^^. In order to obtain second rank 
symmetric tensor from the third rank tensor, which is now the third derivative of (f) like V^Vi/Vp(/), we need 
one more derivative of (j) ^ in the second term of (fTU]) like V^(/)VpVi/Vp0, V (V^V,y Vp(/< -I- V,yVpVp0), or 
V iy fj,^ ufl) + Vi,VpVp0). This kind of terms was not considered in [5:]. 

4. The expression does not contain hu so that the traceless and transverse conditions Cj-* = d'^Cij — for Ctj in 
should be kept as in the Einstein gravity. This condition ensures the absence of the extra modes. 



5. The expression should contain the graviton, which is the traceless and transverse part Cy in ([?7|) of the perturbed 
metric. 



The combination of the terms in ([10)) satisfies the above conditions but the combination in ([10)) is not, however, 
unique. More generally, one can propose the combination like 

Q^p = PJ^Pp" [Pp, - -l-ap^/.VVpV.^^ + {ciSfS^'^ + caP^^P/) V (VpV.Vp0 + V.VpVp0) 

+ (c3<5a^ V + CiP^^Pp'') VP (VpVpV,</) + V,VpV^0) + (csg^. + cePp,) (d'^^d'' ^R^,, + U^R) . (50) 

Here ci, C2, C3, C4, C5, and ce are constants. When we consider the fluctuation from the flat metric in ([4t. 

V (VpV^Vp.?!) + V^VpVp0) - V (VpVpV^0 + Vr.VpV^(t>) - -htt,^^ , 

5^05Vi?p, + UoR ^ Uo {d'd^h^j - dkd'' } . (51) 

Then by using ()17p and ()27p with ([25)) and ([57)) . we find the terms with the coefficients ci, C2, C3, C4, C5, and cq vanish 
and therefore these terms do not affect the propagator of graviton although they may give rise some interaction terms. 
Having in mind the above considerations, one may speculate that proposed theory is multiplicatively-renormalizable 
in gauge- fixed formulation (after spontaneous breaking of Lorentz symmetry) . 

In the above analysis, we have chosen Eq. ([S]) as a gauge condition. Then Eq. (IT7)) follows from the constraint 
equation ([2]). As another gauge condition, instead of (O, we may choose (flT)) . When we write the fluctuation of (/> as 

= y/Wot + 6(p , (52) 

the constraint equation ^ gives 

Then 6(p does not depend of the time coordinate: 



(50 = $(a;) . 



(54) 
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The gauge condition (jl7|) has the residual gauge symmetry for the diffeomorphism with respect to the time coordinate 
t. The residual gauge symmetry is the shift of t by a function independent of t 

t^t + f{x). (55) 

Therefore by using the residual gauge transformation given by 

fix) = -^, (56) 

we obtain ([5]) from ([5^ . 

Thus, we presented gauge-fixed formulation of the theory which leads to only propagating higher derivative graviton. 
Scalar and vector degrees of freedom do not propagate. 

III. SPONTANEOUS LORENTZ SYMMETRY BREAKING 

Due to the constraint equation vector has nontrivial value and therefore the Lorentz symmetry is broken. 
Note that the Lorentz symmetry breaking is spontaneous. Let us compare this situation with that in the usual field 
theory like the Higgs model. The usual U(l) Higgs model, whose potential is given by 

Vm,,s = ~^rcb+^i<P*^f , (57) 

has a global U{1) symmetry, which is the invariance under the transformation 

(/)^e*^°0, (58) 

with a constant real parameter Oq- In (|57p. is a complex scalar field and m and Aq are positive parameters. The 
minimum of the potential is given by 

(59) 



Ao 

Here tp is a constant phase. The value of if can be arbitrary. If one chooses specific value of tp, the value of ip is 
changed under the U{1) transformation ([55]) as ip ^ f + 6, and therefore the ground state is not invariant under the 
U{1) transformation (|58p and the U{1) symmetry breaks spontaneously. One can always choose the real axis of the 
complex (/)-plane to be parallel with the value of (f) in the ground state so that tp — 0. 

In our model, the constraint equation ([2]) shows that the value of the vector (dfj^cj)) is located on the hyperboloid 
defined by 

- xf'x^ = f-x^ = 2Uo . (60) 

The value of the vector {df^cf)) changes on the hyperboloid under the Lorentz transformation. If we choose a value of 
{d^cf)), the Lorentz symmetry is broken spontaneously. After that one can always choose the time axis to be parallel 
to the vector {d^(j)). 

We should also note that the actions pT|) and ((T^ have a shift symmetry 

(f) ^ (P + (f>o . (61) 

Here is a constant. In the flat vacuum background, the actions (jlip and (jl2p are invariant under the time 
translation; 

t^t + to. (62) 

Here to is a constant. Since the solution ([5]) of the constraint equation ([3]) depends on the time coordinate t, the 
solution spontaneously breaks the symmetry under the time translation (|62p . The solution ([5]) also breaks the shift 
symmetry in (pT|) . Note that the diagonal symmetry of the time translation (p^ and the shift symmetry (I5T|) is not 
broken. In fact, if (po in (pT|) is chosen as 

^o = -V^to, (63) 

under the simultaneous transformation, the solution ([5]) is invariant. The diagonal symmetry effectively plays the role 
of the time translation and the flat vacuum solution is effectively invariant under the time translation. 
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IV. VECTOR FIELD FORMULATION 

Let us reformulate the theory in terms of vector field. Note that the actions ([TT|) and (|12p can be rewritten by 
using the vector field A^. For the vector filed A^, if we impose the constraint 



the vector field becomes pure gauge field and can be rewritten as 
Furthermore, by introducing a new scalar field ip and considering the combination 
the combination A'^ is invariant under the gauge transformation 



A., 



Af, + a,,e(a;^) , (p->- ip-e. 



(64) 
(65) 
(66) 
(67) 



Then we rewrite the actions ([Tl]) and ([T2|) . which are invariant under the gauge transformation ([67]) as 



R 
2^ 



2[/n 



1 



-A( -A)JA'^'' + C/o) -B'^'^F^, 



for z = 2n + 2 model (n = 0, 1, 2, • • • ), and 



(68) 



S-. 



2n+3 



|(^^Myi^-VpV,.-2;7oV''Vp)"+'P'^"'^P'^^'^ (^Pp, - J^A^VV^A^^ I 



(69) 



for z = 2n + 3 model (n = 0, 1, 2, • • • ). Here B^'^ is the anti-symmetric tensor field whose variation gives a constraint 
In the actions dlHl) and dSHl), P/^"" is defined by 



-"/^ + 2[/o 



(70) 



instead of ([5]). 

The actions (|68l) and ([M]) are classically equivalent to the actions (ITT|) and p^ . To be sure, let us check the 
propagating modes on the flat background. For this purpose, we decompose P^^ and yl^ as follows. 



Then one gets 



Bu = d j + e^ , Py = diTUj - dj-rrii + eij^d'^k , A,; = 9^ + k ■ 
d'e, = d'm, = d'l, . 



B'^'Ff,^ = 2fdkd'' {-dt(l> + At) - 2m^dkd''lj + 2e^dtli + total derivative terms . 



(71) 



(72) 



Note that in the above expression and therefore even in the total action, there does not appear and k. The variation 
over / gives 



(73) 
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and the variation over gives 

h=0. (74) 

The variation over gives dth = 0, which is consistent with (|74p and does not lead to new constraint. On the other 
hand, the equation given by the variation over Aq can be solved with respect to / and therefore / becomes auxiliary 
field. Similarly, the equation given by the variation over li can be solved with respect to rrii, which also becomes 
auxiliary field. Then from (|7T|) and (|73|) . we find is pure gauge field = dfj,(p. Hence, in the actions (|68p and 
()69|) . the scalar mode (f) and the scalar field ip appear only in the combination of 9^ {(p + if). Therefore the variations 
over (j) and ip give the identical equations corresponding to (j20p and (|2ip. Hence, using the same arguments as in 
Section ini we obtain the same graviton propagator (|47|. 



V. FRW COSMOLOGY 



In this section, we discuss simple accelerating FRW cosmology for higher-derivative gravity under discussion. We 
start with a little bit different but general action: 



S 



R 

2^ 



« |(9'^0-9>VpV, - a^^a^^vvp)" p/p/ (^R^, - ^ap^z-vvpy.^) 



-A ( -ap0a^0 + c/o 



(75) 



In the low energy, one only needs to consider the contribution from the first Einstein-Hilbert term. Since the n = 
term contributes in the next-to-leading order, we now consider the simplified model given by the Einstein-Hilbert 
term and the n — term: 



5* 



d x\/~g 



R 



2k2 



1 



2Uo 



- aoP^^p," ( Pp. - Trfrap^vVpV.^ ) p^A'p/^- ( _ _a^^v''VpV, 



1 



2C/o 



(76) 



It is interesting that even in the FRW background, we can choose the local Lorentz frame and obtain the same graviton 
propagator ([i7|) at short distances, where space-time can be regarded to be fiat. 
In order to consider the cosmology, we assume the FRW-like metric: 



ds' = -e"'^'^dt' + aitf J2 (^^' 



(77) 



i=l,2,3 



Here a{t) is called a scale factor. The variation of the action with respect to b{t) gives the equation corresponding to 
the first FRW equation and the variation with respect to a(t) gives the equation corresponding to the second FRW 
equation. Since we have 



^ it — ^ ti — ^ tt — ^ jk — ^ ■) 



(78) 



and 



Rt 



3{H + H^ -bH 



Ri 



H + ZH^ - hH a^e-^^Sii , R,t = Ru = , 



R^6[H + 2H^ -bH]c 



.-2b 



(79) 
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the action (l76l) has the fohowing form: 



(80) 



We should note dt4> = e^y/2UQ and di(j) = 0. It is remarkable that this action is very similar to the one of well-known 
gravity. By the variation of the action with respect to b, one obtains an equation corresponding to the first FRW 
equation: 

+ 81aoH^ = Pmattcr • (81) 

Here /Omattcr is the energy-density of the matter which was not explicitly written before. One puts 6 = after the 
variation over h. By the variation over a, one also obtains the equation corresponding to the second FRW equation, 
which can be evaluated by using the first FRW equation and matter conservation law. We should note that the equation 
([STjl is the first order differential equation with respect to the scale factor a{t), which should be distinguished from 
the usual higher derivative gravity like F{R) gravity, where the equation corresponding to the first FRW equation 
is the third order differential equation with respect to a{t). Note that H = 0, which expresses the flat solution, is a 
trivial solution of (ISTI) when there is no matter. We also note that there is a de Sitter solution, where H is a. constant, 
given by 

H' = -:^^, (82) 



which may express the inflation in the early universe. The existence of the de Sitter solution in (|82p requires ao < 0, 
which does not always conflict with the condition to avoid the tachyon in (|48|) . Now we only kept the term with n = 
in the action (j75p . just for simplicity. If we include the higher term with n > 0, the de Sitter solution might become 
unstable and there can occur the instable inflationary solution. 



VI. DISCUSSION 



In summary, we formulated covariant higher derivative gravity with Lagrange multiplier constraint and scalar 
projectors. It is demonstrated that such theory admits flat space solution. Its gauge-flxing formulation is fully 
developed. The study of spectrum shows that the only propagating mode is (higher derivative) graviton, while scalar 
and vector modes do not propagate. Eventually, scalar and vector modes correspond to composite states at any 
perturbative level. 

Furthermore, we show that Lorentz symmetry breaking in the theory under discussion is spontaneous. The equiv- 
alent formulation in terms of vector flelds is developed. The preliminary study of FRW cosmology indicates to the 
possibility of inflationary universe solution. It is interesting that flrst FRW equation in the theory turns out to be 
the first order differential equation which is quite unusual for higher derivative gravity which normally leads to third 
order differential equation with respect to scale factor. This may indicate to presence of some hidden symmetry in 
the higher derivative gravity under consideration. 
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